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Abstract — Measuring of thin film properties is difficult when 
compared to bulk materials. A straight, horizontal cantilever 
beam exposed to free vibrations will vibrate at its characteristics 
or natural frequencies. In this research paper the vibrations of a 
thin film cantilever beams will be studied and the equations of 
vibrations will be derived. The frequency equations are solved 
using Matlab® to show the output frequencies and the mode 
shapes related to each frequency 

Index Terms — cantilever, the mode shapes 

I. Introduction 

A cantilever beam is one of the most fundamental structural 
and machine components used in many different applications 
for decades. Cantilever beams are generally beams with one 
end fixed and the other end free. The length has a much larger 
dimension when compared with the width and depth. In 
addition, cantilever beams maybe straight or curved, with 
rectangular or circular cross sections. Figure 1.1 shows a 
cantilever beam with a rectangular cross section. 



Figure 1.1: Rectangular cantilever beam 


Cantilever beam design and shape depends on the 
application. Its size, material, and weight are different from 
one application to another. For example, one of the most 
common applications of a cantilever beam can be shown in 
figure (1.2) below. The “fixed wing” in meters is designed as 
a beam for some preliminary analysis to help lift the plane 
and make it fly [1] 

Wings of a Stationary Aircraft is a Cantilever Beam 



Figure (1.2) Source : Student Online Laboratory through 
Virtual Experimentation 

In Micro Electrical Mechanical Systems (MEMS), micro 
cantilever beams are used in radio frequency filters and 
resonator. Beams are exposed to different dynamic loads in 
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different applications. When a cantilever beam is exposed to 
a dynamic load, the load will excite the beam to vibrate at its 
characteristic, or natural, frequencies. Studying vibrations of 
cantilever beams are very important. It helps to determine the 
durability concerns (by analyzing dynamic stresses) and 
noise .This information used to reduce the discomfort and 
excessive stresses in different applications in which beams 
are essential components. 

Cantilever beams can be more than one layer, in many 
applications beams are coated with one or more layers. This 
is possible especially in MEMS applications. In these 
applications, we may need three layers. One of the MEMS 
examples is a beam coated with a piezoelectric material. This 
beam is usually called a “sandwich beam”. In this beam, two 
piezoelectric layers are used to coat the beam from both 
sides. As a result, exposing this beam to vibrate will lead to 
excite the piezoelectric materials to generate a voltage 
difference between the two piezoelectric layers. This voltage 
can be used in many applications to measure stress, strain, 
and get many useful outputs. Another example is a two layer 
cantilever beam, in which a beam is coated with a different 
material. For instance, an aluminum cantilever beam is 
coated with thin film zinc oxide to form a two layers 
cantilever beam. The coated Aluminum beam has different 
characteristics and frequency responses than single 
Aluminum beam. Frequency values will be changed after the 
beam has been coated. These changes can be used to measure 
different characteristics. For example, the modulus of 
elasticity for thin films cannot be measured using 
conventional methods, because measuring modulus of 
elasticity for zinc-oxide is very hard and almost impossible. 
This is due to the thickness of the film and its sensitive 
characteristics. Thus, using vibration analysis will help in 
measuring the modulus of elasticity for thin film materials 
[ 2 ]. 

As seen above, many applications exist for multi-layered 
beams. In this research paper, we will study the frequency 
behavior for a one- and a two-layer cantilever beam. The 
equations of motion, frequency and characteristic equations 
will be derived. Numerical results will be calculated to see 
the frequency differences between one and two layered 
beams. Different materials will be taken to compare the 
results. In addition, results and future work will be discussed 
at the end of this paper. 

II. LITERATURE REVIEW 

In this research our concerns are to calculate and explain the 
free vibrations of two-layer cantilever beams. However, to 
make it more obvious we will explain the vibration of a one 
layer cantilever beam before we move to a two -layer 
cantilever beam [2] . 

In this way, we can see the difference in equations and 
characteristics when we move from one layer to two layer 
cantilever beams. Our calculations in both cases will depend 
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on free and un-damped vibration. Thus, the beam isn’t 
exposed to any external force. According to Gorman (1975) 
“In free vibration a beam undergoes oscillatory motion while 
free of any external forces, whereas in forced vibration the 
beam responds to a system of time varying external forces”. 
([2] -[4]). Free vibration of cantilever beams can happen in an 
infinite number of mode shapes, each mode has a discrete 
frequency. 



The first frequency which is the lowest one is associated with 
the first mode; the second frequency is associated with the 
second mode and so on. However, higher frequencies - third 
and above - are less significant. This is because they are 
difficult to excite and the number of points on the beam 
having zero displacement increase directly with the mode 
number. To get appreciable amplitude for the higher modes, 
much more energy is required in this case [5]. We will take a 
bulk material coated (covered) with thin film in our analysis 
even though the analysis will be valid for any kind of beams. 
A thin film is a layer of material ranging from fractions of a 
nanometer (monolayer) to several micrometers in thickness, 
and usually the act of applying a thin film to a surface is 
called thin-film deposition. For example, an Aluminum 
cantilever beam can be coated with Zirconate Titanate (PZT) 
thin film to form a two layers cantilever beam. Dealing with 
thin film materials needs more attention because thin films 
characteristics are changed rapidly, with non-linearity, and 
sensitivity. One important thing about dealing with double 
layer cantilever beam that the thin film beam stiffness should 
be different from the main beam (bulk material), so the 
frequency shift can be noticed ([5] -[7]). 

Two main methods are used to obtain the solutions for a free 
vibration cantilever beam. The first one depends on solving 
differential equations for the equilibrium between inertia 
forces and elastic restoring forces subject to boundary 
conditions. The second method called the energy method, 
this method depends on the fact that the sum of the potential 
energy and the kinetic energy is always constant. [2]. 

In the next chapter we will use the differential equation 
method to derive the equations for free vibration of a 
cantilever beam. These equations will be the main reference 
for the next chapters in which it will be used to derive the 
equations for two layers beam, and obtain the frequencies 
and mode shapes. 

III. CANTILEVER BEAM FREE VIBRATION THEORY 

A. Equation of Motion 

We assume a cantilever beam with length at least 20 times the 
average depth, and the beam vibrates transversely in the 
z-direction as shown in figure 2.1 below. It makes the 
problem mathematically one dimensional with no torsional 
vibrations. We take a small differential beam element of 

length dl from the beam as shown in figure 2.2. By 
developing the free body diagram for this element, we can 
see that there is a shear force and a bending moment in both 
sides. These quantities vary along the beam with time (t). In 
addition, in free vibration case there is no force, p=zero 
([l]-[3]). 



M+dM 


V+dV 


Figure 2.2: Differential 


Figure 2.1: Cantilever beam 
beam element 

The summation of forces around the center line is: 

d 2 w 


X F >’ = 


m 


at 



Where p is the mass per unit volume and A is the cross 
sectional area, the beam is assumed to be under small 

dw 

vibrating amplitude; thus, the slope of the centerline 


is 


dx 


d d 

everywhere small and — — — as well. Cancelling terms 

dl dx 

and dividing by dx (2.1) becomes: 



— > ( 2 ) 

Take the moments about the center: 

0 

dM 

—M + (M + dx)-Vdx = 0 ->( 3 ) 

dx 

Cancelling terms and divided by dx , equation (2.3) 
becomes: 


-►( 4 ) 


Back to Euler-Bernoulli equation for linear elastic materials: 




— » ( 5 ) 


Where 


d 2 w 


is the small slope curvature, E is the Young’s 


modulus of elasticity, and I the second moment of inertia 
area. Substituting (5) into (4) gives: 


— > ( 6 ) 



Now substituting (6) into (2) gives: 
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) + pA 





->(7) 


This is the equation of motion for homogeneous, constant 
cross section, and free vibration cantilever beam. Next, we 
will solve this equation for free vibration cantilever beam. 


IV. Solution of the Differential equation 


Using the method of separation of variables we assume a 
solution of the form: 


v*Z,r) = X(Z).G(f)->( 8) 


X: independent on time , Q: independent of position 
We rearrange equation (7) into two parts, one depends on 
position and the other depends on time: 



)) 

dl 4 


d 2 (X(l).Q(t )) 
= -m 

dl 2 


Divide two sides by m.X(l).Q(t ) we get: 

— > (9) 


2 

Each side is constant and equal to CO which is the natural 
frequency of the cantilever beam with real quantity, take the 
left side of the equation and rearrange it to become: 

d 4 X(/) rmf 

dl 4 El 

P 


/? 4 = The corresponding non-dimensional frequency for 
infinite number of values. To determine /? we will find the 

general solution for the differential equation (8). 

The general solution of equation (2.8) is given by: 





EI d 4 X(l) _ i_ a 2 g(p 

mX (/) dl 4 ~ Q(t )' dt 2 

\ „ j \ „ > 


-►( 11 ) 


Because we are interested in finding pi , we will calculate the constants Ci C2, C3, and C4 for equation (11) only. To find 

these constants and solve the equation, boundary conditions and initial conditions should be determined. However, the above 
equations are valid for any beam with arbitrary boundary conditions. In our case we will solve these equations using the 
boundary conditions for a cantilever beam [8] 



V. Boundary conditions for a cantilever beam in free vibration 


In this section we will derive the solution for the cantilever beam shown in figure 2.3 below. This 
beam has a length (L), a thickness (t), and a width (w). 

The solution here is for the one layer cantilever beam. We will use this solution to derive the solution for two-layer cantilever 
beams. Applying the following boundary conditions for this beam: 
a- At the fixed end the deflection (z) and the slope are zeros. 

w(0,f) = 0-^(0,*) = 0->(13) 

cl 


b- At the free end the moment(M) and the sheer force(V) are zeros 


d 2 w d 3 w 

— — (L, 0 = 0^ — — (L, 0 = 0^ (14) 

dl 2 dl 3 



Figure (2.3): One layer cantilever beam 

Implement equations (13) and (14) into equation (8), we get X ( l ): 
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We will use the boundary conditions (15) to solve equation (1 1) by differentiating it three times. 

->(16) 



d~X 

~df 


= / 3 2 (-C , sin (/?/)- C 2 cos (/?/) + C 3 sinh (/?/) + C 4 cosh (pi)) — >• (17) 


, <fx_ 

~df 


= J3 3 (- C,Cos(j3l ) + C 2 sin (/?/) + C 3 cosh (01) + C, sinh {pi)) -» (1 8) 


Substitute the boundary conditions in equations (11) and (16) we get : 

X (0, t) = 0 = Cj siX(0) + C 2 cos(0) + C 3 sinXiO) + C 4 cosh(0) 


C = —C 


(19) 


dX 

~df 


i l=0 = 0 = p(C\ cos(0) - C 2 siMO) + C 3 cosh(0) + C 4 sinX(0)) 


C, = -C 3 (20) 


Now substitute equations (19) and (20) with boundary conditions into equations (17) and (18) to get: 



->( 21 ) 


—>( 22 ) 


These equations can be written in a matrix form: 


^ sin(/?L) + sinh(/?L) cos(/?L) + cosh(/?L) V 


v 


cos(/?L) + cosh(/?L) - sin (J 3 L) + sinh (/?L) 


J 


i 


\ C 2J 


^0^ 

vOy 


To solve this matrix we will take the non-trivial solution at the resonant frequency in which the determinant of the matrix must 
equal to zero; thus: 

sinh 2 (fiL) - (s/n 2 (fiL) + c0s 2 ( /3L )) - 2 cos (J3L) cosh (J3L) - co^h 2 (J3L) = 0 

1 + cos (J3L) cosh (J3L) = 0 — » (23) 


(3 > 0 , p ^ 0 (No vibration when the beam at rest) p ^ negative (Negative values lead to redundancy) 

This is a non-linear equation with infinite number of roots (i.e. continuous system). Solving this non-linear equation (** See 
Appendix-A) we get: 

PL — 1 .87 5, 4.6904, 7.548, 

Back to equation (2.10) and rearrange it for CO : 

( BLf [W 

® = ~r— . — t ->(24) 

| L \pA\ 

Where: E = Young modulus of elasticity ( Pa ) 

1 

/ = —bt = Moment of inertia- Area for a rectangular cross-section ( m ) 

p = Mass density ( Kg / m ) 

A-tb= Area of the cross section ( m 2 ) 

L,b,t = Length, width, and thickness respectively 

This equation computes the frequencies for a one layer cantilever beam. Infinite number of frequencies can be calculated using 
this equation. In the next section we will derive the equation for a two-layer cantilever beam [10]. 


VI. Vibration of a two layer cantilever beam 

We discussed the vibration for one layer cantilever beam, but 
many recent applications use more than one layer in which 
the beam is coated to become a composite material. 


We will get a new material that is different from the old 
material. This is very helpful in many applications like 
automotive engineering, MEMS and NEMS, and biomedical 
equipment. The use of composite materials can reduce 
weight, increase efficiency and improve output response. As 
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a result, instead of using Steel for a cantilever beam, we can 
use two other materials (composite) and get the desired 
output. In our case get the same frequency response. 

Thin film materials usually in micros-meters or even 
nano-meters are using to coat other materials, in this case we 
will take a beam coated with a thin film to form a new two 
layer beam as shown in figure (2.4) below [9]. The new beam 
will have new characteristics and new specifications. This 
two layer beam has two materials. Each material has different 
density, and Young modulus of elasticity. In addition, inertia 
and thickness differ from each layer to the other. Both layers 
still have the same length and width, and both of them will 
vibrate at the same frequency. 


A 1 v 12 v 21 v 23 v 32 v 31 v 13 2 v 21 v 32 v 13 


V 12 , V 21 , V 23 , V 32 , v 31 , V 13 = The poisson’s ratios 

Where I n is the modified mass per unit length, is the 

material Elasticity, and N is the number of layers. These are 
the general equations for any number of layers, in our case for 
a two-layer beam, we will substitute N=2, so equations (25) 
& (26) become: 

I, 


k = 1 



-+ (27) 



Figure (2.4): Thin film cantilever beam 
Both materials will vibrate at the same frequency, so the 
frequency will be shifted to a new value (i.e. modified). This 
frequency shift will be very small, because the thin film is in 
micro-scale. We will take the case of a two layer beam and 
derive the necessary equations. More layers can be added 
easily and implemented to the equations. 

Lets take the beam in figure (2.4), where (f) stands for the 
thin film and (b) stands for the coated beam (bulk), the beam 
has: 

3 

p h = Density of the bulk material in ( Kg / m ) 

Pj = Density of the thin film material in ( Kg / m ) 


E b = Young Modulus of Elasticity for the bulk material 

(Pa) 

E j = Young Modulus of Elasticity for the thin film 

(Pa) 

tj,t b = Thin film thickness and Bulk material thickness 
respectively (m) 

Back to frequency equation (2.24), the modulus, inertia, and 
density will become: 

— » (25) (Qatu &Leissa, 2011, 

p375) 

— » (26) (Qatu &Leissa, 




2011, p375) 



( 1 v 23 v 32 ^ 


V 


A 


J 


(Qatu 


Figure (2.5): Two layers cantilever beam thickness 
The new modulus of elasticity will be: 

D„=l>e;, 

k = 1 






) 


Because we are dealing with one direction only and for more 
simplicity, we consider ADI and V 23 , V 32 =0 .Thus, the 
equation becomes: 

-+(28) 

Where E is the modulus of elasticity of each material, t is the 

layer thickness, and b is the beam width. 

We can implement these equations now in the frequency 
equation (2.24) we obtained earlier 



co jnf!K 

L- V /, 



This is the frequency response for a two layer cantilever 
beam. We can add more layers by following the previous 
steps for two layers. 


VII. Properties of the mode shapes 


&Leissa, 2011, p367) 

After we define the equations for the frequency for both cases (one layer and two layer beams, 


we need to know the mode shapes of the frequencies. 
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There is infinite number of non-dimensional frequencies and eigen values for each problem. Back to equations (2.1 1), we can 
rearrange it for the eigen vectors to determine the corresponding eigen functions. (Leissa & Qatu, 2011, p.l 10), the equation 
will become: 



X (/) = Cj sin (/?/) + C 2 cos(/?/) + C 3 sinh(/?/) + C 4 cosh(^/) 

Using equations (2.19), (2.20), (2.21) and (2.22) for a cantilever beam and combined them we get: 


X(l) = 

cosh(y^/)-cos(^/) 
sinh(/?/) - sin(/?/) 

. [sinh (/?/) - sin (/?/)] + cosh (/?/) - cos {fil) 


(29) 


->(30) Where 


t = _ c°shW-cosW = o 734> i 01 8 , 0 . 9 99, 

sinh(/?/)-sin(/?/) 


So equation (30) becomes: 



X (/) = the characteristic function for each frequency 
Back to equation of motion (8): 

v*/,f) = x(/).G(0-K 8) 

This equation satisfies the boundary conditions in (13) and (14), which is the solution for each frequency; however, the general 
solution of the equation of motion is the sum of the characteristic vibrations (Volterra,1965, p.3.18) , thus (8) becomes : 

OO 

<7,0 = Y^ X n^)\- A n cos( w n t) + B n sm(rv,/)] 


W 


n = 1 


W 


(/, t) = Yj\- C0Sh( ^! ) —^7 .[si nh( /?/) - sin (/?/)] + cosh (/?/) - cos(fil ) [4 cos(w/) + B n sin(iv„/)l 


n = 1 


sinh(/?/)-sin(/?/) 


*The above equation (32) 


is the general solution for any number of modes (n) 


To determine the coefficients A n and B ri , we will use the initial conditions 


n 



Both parameters are constants and depend on these conditions 

w(l,t) = (p 

Thus, Qw 


dl 




—> where 0 <1 < L 


It follows from equation (32) above: 


oo 


Y J A n X,(D = <p(l) 


n = 1 

oo 


o=ao 


n = 1 


- > ( 33 ) 


These equations are analogous to the Fourier series expansions. Thus, both equations are determining the series expansions of 
the functions in terms of the characteristic functions X n (Volterra & Zachmanoglou, 1965, p.318). 

Finally, to find the constants A n and B n we make use of the orthogonality for the characteristic functions; thus, by 

integrating these equations with respect to (/) fromO — > L , and multiplying by X (/) we get the formulas for A and B : 
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— > ( 34 ) (Volterra & Zachmanoglou, 1965, p.913) 


In the last equation we can eliminate B n , because it depends on the initial velocity, and in this case the beam starts at rest. 

Thus, B u =0. So for our calculations we will evaluate A n , because A n depends on the initial position at t=0 only. 

In the next chapter we will evaluate the displacement caused by each mode, and the total initial displacement for different 
cantilever beams, as long as we will draw the mode shapes for each case using MATLAB program. 


VIII. EXAMPLES AND RESULTS 

Material Properties and Selection 


We will introduce first some common materials properties. 
Determining the suitable material is very important in any 
cantilever beam design, because the material properties play 
a significant role in the output values and the whole beam 
vibration. Thin films are sensitive to circumstances and 
different external disturbances. Their properties are changed 
rapidly with instability. For instance, the density for a thin 
film material can be easily affected by temperature, and that 
results in changing the outputs and the related frequencies as 
well ([10]-[12]). 

In our calculations of beam vibration, we need to determine 
three main characteristics: beam dimensions (length, width, 
and thickness), Young Modulus of Elasticity, and material 
density. All of them are dependent on the application in 
which the beam will be used. For instance, in our experiment 
we will use a thin film cantilever beam, a beam coated with a 
thin layer; the dimensions will vary from micros to 
centimeters, also the density and the Elasticity should be 
chosen carefully. 

In the next table (3.1) we summarize some of the materials 
specifications. We concentrate on the density and the 
elasticity which is important for our calculations. 

Table (3.1): Materials properties 


Material 

Density 

(kg/m A 3 

) 

Young 

Modulus 

(GPa) 

Reference 

Aluminum (AL) 

2700 

70 

CRC Materials 
Science and 
Engineering 
Handbook, p.46 

Silicon (Si) 

2330 

190 

J.mater.Res,Vol.l2, 
No.l,Janl997, p.59 

Copper (Cu) 

8960 

128 

Mechanics of 
Materials 23(1996), 
p.314 


Lead Zirconate 
Titanate (PZT) 

7750 

10-40 

IEEE Micro Electro 
Mechanical 
Systems Workshop, 
Jan-Feb 

1 99 l,Nara, Japan, 

p.118 

Stainless 

steel(Bulk) 

7900 

200 

Proceedings of 
IEEE,Vol 

70, No. 5, May 1982, 
p.421 

Stainless steel 
(Thin film) 

7900 

240 

Thin solid films 
290-291(1996), 
p.363 

Platinum 

21440 

170 

IEEE, Micro 
Electro Mechanical 
Systems Workshop, 
Feb 1990, Napa 
Vally, California, 
p.174 

Silicon Nitride 
(SixNy), thin 
film 

2843 

-2887 

85-105 

Sensors and 
Actuators 
A, 35(1992), 
p.153-159 

Titanium (Ti) 
(Bulk) 

4510 

115 

CRC Materials 
Science and 
Engineering 
Handbook, p.47 

Titanium (Ti) 
(Thin filn) 

4510 

110 

Thin Solid 
Films, 270(1995), 
p.263 


It is very important to notice the difference in properties for 
the same materials. This table shows the common materials 
used in laboratories and experiments. For the best results we 
will choose the material with stable elasticity and 
density.These characteristics in the table are based on 
experiments, computer analysis, and other technologies; 
thus, there are differences in the same material. Also, as we 
note that materials may have two categories, one for bulk and 
the other for thin film. For instance, Stainless steel has a 
Young Modulus of 200 GPa for the bulk material, and 240 
GPa for the thin film. 

IX. Calculations and Results 

In this section we will calculate the frequencies and mode 
shapes for one and two layers beam with different materials 
selection. We will use MATLAB to implement the equations, 
calculate the frequencies, and plot the mode shapes [12]. 
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The program will calculate the frequencies in Hertz, and can 
plot up to six mode shapes. The first modes have larger 
values and the higher vibration modes are not shown. 

• Aluminum cantilever beam (one layer) 

For Aluminum we consider the following characteristics: 

L= 0.12 m b=0.03 m 1=0. 00 2 m 

E=70 GPa p =2700 Kg / m 3 

1. Natural Frequencies (Hz) 

Mode shape ( 1 ) (Hz) : 113 .420 
Mode shape (2) (Hz): 710.792 
Mode shape (3) (Hz): 1990.239 
Mode shape (4) (Hz): 3900.074 
Mode shape (5) (Hz): 6447.102 
Mode shape (6) (Hz): 9630.854 

2. Mode shapes 


2. Mode shapes 


Mode shapes of the Cantilever beam 



Figure (3.2): Steel beam mode shapes 

• Aluminum coated with Zinc oxide thin film (two layers 
beam) 

For this beam we consider the following characteristics: 


Mode shapes of the Cantilever beam 



L h — L f — 0.12m b b — b f — 0.03m 
t b = 0.002m tj — 0.00002m 
E b = 70 GPa E f = 1 1 lGpa 
p b = 2700 Kg / m 3 p f = 5680 Kg / m 3 

1. Natural Frequencies (Hz) 

Mode shape (1) (Hz): 1 14.933 
Mode shape (2) (Hz): 720.272 
Mode shape (3) (Hz): 2016.782 
Mode shape (4) (Hz): 3952.088 
Mode shape (5) (Hz): 6533.085 
Mode shape (6) (Hz): 9759.297 

Notice the slight difference after the Aluminum coated, 
because the thickness of the thin film is only 20 micrometers. 

2. Mode shapes 


Figure (3.1): Mode shapes for Aluminum beam 
Figure (3.1) below shows the displacement caused by each 
mode at t=0 as long as the total displacement for the beam. 
From these displacements it is obvious that modes higher 
than the third mode have a little effect on the total 
displacement. 

• Steel cantilever beam (one layer) 

For Steel we consider the following characteristics: 

L= 0. 1 2 m b=0.003 m 1=0.002 m 

E=210 GPa p =7850 Kg / lit 

1. Natural Frequencies (Hz) 

Mode shape (1) (Hz): 116.044 
Mode shape (2) (Hz): 727.235 
Mode shape (3) (Hz): 2036.279 
Mode shape (4) (Hz): 3990.295 
Mode shape (5) (Hz): 6596.244 
Mode shape (6) (Hz): 9853.646 

The frequencies depend on the Young Modulus and the 
density for the material. 


Figure (3.3) below shows the mode shapes for a coated beam 
(two layers). It is so difficult to notice the changes. In 
addition the mode shapes amplitude is the same, in our work 
we don’t care about the amplitude, because we don’t have 
any damping system to change it; thus, it will be the same in 
all cases. 


Mode shapes of the coated Cantilever beam 



Figure (3.3): Aluminum coated with Zinc oxide mode shapes 
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• Steel coated with Silicon nitride (two layers beam) 

For this beam we consider the following characteristics: 

L b —L f — 0.12m b h — b f — 0.03m 
t b — 0.002m tj — 0.00002m 
E b = 210GPa E f = 190Gpa 
p b = 7850 Kg / m 3 p f = 2330 Kg / m 3 

1 . Natural frequencies (Hz) 

Mode shape ( 1 ) (Hz) : 117.451 
Mode shape (2) (Hz): 736.054 
Mode shape (3) (Hz): 2060.973 
Mode shape (4) (Hz): 4038.684 
Mode shape (5) (Hz): 6676.234 
Mode shape (6) (Hz): 9973.138 

2. Mode shapes 


Mode shapes of the coated Cantilever beam 



The last examples show the behavior for different cantilever 
beams according to the frequencies and mode shapes; the 
MATLAB which is shown in (*Appendix-B) computes the 
frequencies and mode shapes for any kind of cantilever beam 
from bulk materials to thin films. 


X. SUMMARY AND FUTURE WORK 

This research is concerned with the frequencies and mode 
shapes of a thin film coated cantilever beam exposed to free 
vibrations. Equations of frequency and mode shapes are 
derived and solved for different cantilever beams with 
different materials. For instance, it explained different cases 
for one and two layer beam, with different materials for each 
case. We used Mat-Lab® to find the frequencies and the 
mode shapes related to them, and compared the results for 
each case. This program is valid for any kind of beams with 
one or two layers and any kind of materials, the user select 
the program for one or two layer beam, and then enter the 
beam dimensions and properties, in a second the program 
calculate the frequencies and the attached mode shapes for 
each frequency. 

Studying the vibration for composite materials (i.e. two and 
more different materials) is very important for a lot of 


industries including manufacturing companies, and 
automotive. The cantilever beam is the main component used 
in many applications and industrial products. Automotive 
companies are always search for lightweight materials in 
their cars; lightweight materials and thin films are highly 
recommended for automotive industries because: 

1 . It has less exhaust and it prevents high pollution 

2. Cost reduction 

3. Recycling and environmental friendly 

4. Maximize fuel efficiency 

5. Reducing vehicle weight 

Automotive industries are searching for lighter, stronger, 
safer, and more environmental friendly materials than current 
ones. However, studying the frequency responses for the 
materials is very important for automotive industries as long 
as the material itself. The automotive companies are 
concerned with reducing vibrations on the vehicle. 
Vibrations affect the vehicle durability, operation, and 
stability. Thus, knowing the vibration will help to make 
improvements and increase the vehicle durability. Beams 
generally and cantilever beams specially are one of the main 
components of the vehicles; thus, studying and analyzing the 
vibrations of these components are very important in the long 
term as long as the short term. 

A future work will be on measuring the frequencies for real 
beams in the laboratory and comparing the theoretical results 
with the practical ones. In addition, it is possible to find 
material properties that it’s hard to find in the basic ways, and 
see how the frequency changes with external disturbances 
and effects. 
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• APPENDIX-A: SOLVING NON-LINEAR EQUATION 

Solve the non linear equation 1 + COS(/?JL) COSh(/?Z^) = 0 

Let X- [51 

cos(x) cosh(x) = -1 

From the identity 

cosh(x + iy) = cosh(x) cos(v) + i sinh(x) sin(v) 

- 1 + i'O = cosh(x) cos( y) + i sinh(x) s/n( y) 

-l = cosh(x(l + /)) 

Let Z = x(l + l) 

From the identity 


cosh(z)= -1 



e 2 + e~ z = -2 


(e z ) 2 + 1 = -2e' 


(e z ) 2 +2e z +1 = 0 ^ (e z + IY =0 -+■ e z =-l 


Use the identity 


e lK - cos (n) + i sin(/r) 


— > e z = e in — > z = in + 2nn 


So, x(l + i) = m + 2 nn —> x - 


7ri + 2mr 
1 + / 


n + 2 nn . n — 2rur 

Take the conjugate — ^ X — h l 


Take the real part — ) 



, for n=0, 1, 2, 3. 


• APPENDIX-B :MAT LAB CODES 


I. One Layer Cantilever Beam 


f unctionCbeam (~) 

% Cbeam.m Cantilever Beam calculations 
% HELP: This script computes mode shapes and 

corresponding natural 


o 

o 

clearall ; 
clc; 

closeall ; 
display 

^ I 'k'k~k'k~k-k~k-k~k-k~k'k~k-k~k-k~k-k~k-k~k'k~k'k~k-k~k'k~k-k~k'k~k'k~k-k~k-k~k'k~k-k~k 
-k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k • ^ 

display 

display ('Mohammad Abu-Shams ' ) 
display ('School of Engineering and 
Technology ' ) 
display 

display 

-k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k > ^ 

display ('This program is ONLY for cantilever 
beams with a uniform rectangular cross 
section . ' ) 

display ('This program computes mode shapes and 
crossponding natural') 

display ( ' frequencies up to 6 modes for a 
rectangular cantilever beam.') 
display ('The cross section should be 
rectangular, enter 1 if this is correct, other 
wise the program will be closed') 
display ('To see a quick example , enter 2') 
CS=input ( ' Enter your choice: '); 

9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 919 - 9 - 9 - 

ooooooooooooooooooooooooooooooooooooooooooooo 

9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 

ooooooooooooooooooooooooo 

9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 

ooooooooooooooooooooooooooooooooooooooooooooo 

9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 

ooooooooooooooooooooooooo 

if CS==1 

W=input ( 'Enter Width of the beam in [m] : ' ) ; 
Th=input ( 'Enter Thickness of the beam in [m] : ' ) ; 
L=input ('Enter Length for the beam in [m] : 

' ) ; 

Ix= (1/12) *W*Th A 3; 

A=W*Th; 

disp ('Material properties of the beam') 
disp ('Do you know your beam material 
properties, i.eYoung"s modulus and density ?') 

YA=input ( 'Enter 1, if you do; enter 0, if you 
do not : ' ) ; 

if YA==1 

E=input ('Enter Young"s modulus in [Pa] : 

' ) ; 

Ro=input (' Enter material density in 
[ kg/m A 3 ] : ' ) ; 

else 

display ('Steel: E=2.1e+ll [Pa]; Ro=7850 
[Kg/m A 3 ] ' ) 

display ('Copper: E=1.2e+ll [Pa]; Ro=8933 
[Kg/m A 3 ] ' ) 

display ('Aluminum: E=0.69e+ll [Pa]; Ro=2700 

[Kg/m A 3 ] ' ) 

display ('Zinc-oxide: E=l.lle+ll [Pa]; Ro=5680 
[Kg/m A 3 ] ' ) 

display ('PZT: E=0 . 63-1 . le+11 [Pa]; Ro=7750 
[Kg/m A 3] ') 

E=input ('Enter Young"s modulus in 
[Pa]: '); 

Ro=input (' Enter material density in 
[ kg/m A 3 ] : ' ) ; 

end 

elseif CS==2 

display (' Quick Example') 

display (' This is a rectangular steel beam ') 
display (' Length=0 . 45 m; Width=0.04 m; 
Thickness=0 . 003 m; ' ) 

L= .45; W= .04; Th=.003; 

A=W*Th; 

Ix= (1/12) *W*Th A 3 ; 

E=2 . l*lell ; Ro=7.85*le3; 
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else 



holdof f; 



F=warndlg ( 'WARNING: INVALID 

ENTRY 


elseif MMS==3 



Warning ! ! ' ) ; 



plot (xl , Xnx ( 1 , 

:) , 

' b- ' ) ; hold on 

waitf or ( F) 



plot (xl , Xnx ( 2 , 

:) , 

' r- ' ) 

display (' Type in:>>Cbeam') 



plot (xl , Xnx ( 3 , 

:) , 

' m- ' ) ; grid 

pause ( 3 ) 



title ('Mode shapes 

of the Cantilever beam') 

return 



legend ('Mode #1 , 

'Mode #2 ' , 'Mode #3 ' , 0) 

end 

o o o o o o 

9-9-9-9-9-9-9-9- 

xlabel ( ' x/L ' ) ; 
holdof f; 
elseif MMS==4 

ylabel('Mode shape X n(x) ') 

ooooooooooooooooooooooooooooooo 

O \ /r _ ^ / / OOOOOOOOOOOOOOOOOOOOO 

o/ / rlULlco \ \ ooooooooooooooooooooo 

o o o o o o 

o 

o 

oooooooo 



display ( 'How many modes and mode shapes would 

plot (xl , Xnx ( 1 , 

) , 

' b- ' ) ; hold on 

you like to evaluate ? ' ) 



plot (xl , Xnx ( 2 , 

) , 

' r- ' ) 

HMMS=input (' Enter the number of 

modes 

and mode 

plot (xl , Xnx ( 3 , 

) , 

' m- ' ) 

shapes to computed: '); 



plot (xl , Xnx ( 4 , 

) , 

' c- ' ) ; grid 

if HMMS>=7 



title ('Mode shapes 

of the Cantilever beam') 

disp ( ' ' ) 



legend ('Mode #1 , 

'Mode #2', 'Mode #3', 'Mode 

warning (' WARNING : Up to 6 mode 

shapes 

(plots ) 

#4', 0) 



are displayed via this script') 

• 

r 


xlabel ( ' x/L ' ) ; 

ylabel('Mode shape X n(x) ') 

disp ( ' ' ) 



holdof f; 



end 



elseif MMS==5 



Nm=3*HMMS; 



plot (xl , Xnx ( 1 , 

) , 

' b- ' ) ; hold on 

j j=i; 



plot (xl , Xnx ( 2 , 

) , 

' r- ' ) 

while j j <=Nm; 



plot (xl , Xnx ( 3 , 

) , 

' m- ' ) 

betaNL (jj)=fzero(@ (betaNL) cosh (betaNL) *cos (be 

plot (xl , Xnx ( 4 , 

) , 

’g-' ) 

taNL) +1, [ j j j j+3] ) ; 



plot (xl , Xnx ( 5 , 

) , 

' k- ' ) 

j j = j j +3 ; 



grid 



end 



title ('Mode shapes 

of the Cantilever beam') 




legend ('Mode #1 , 

'Mode #2', 'Mode #3', 'Mode 

index= (betaNL~=0 ) ; 



#4 ' , 'Mode #5 ' 

0) 


betaNLall= (betaNL (index) ) ' ; 



xlabel ( ' x/L ' ) ; 

ylabel('Mode shape X n(x) ') 

%fprintf ( ' betaNL value is %2 . 3f \n ' , betaNLall ) ; 

holdof f 



betaN= (betaNLall/L) '; 



elseif MMS>=6 



k=l ; 



plot (xl , Xnx ( 1 , 

) , 

' b- ' ) ; hold on 

wn=ones (1, length (betaN) ) ; 



plot (xl , Xnx ( 2 , 

) , 

' r- ' ) 

fn=ones (1, length (wn) ) ; 



plot (xl , Xnx ( 3 , 

) , 

' m- ' ) 

while k<=length (betaN) ; 



plot (xl , Xnx ( 4 , 

) , 

’g-' ) 

wn (k) =betaN (k) A 2*sqrt ( (E*Ix) / (Ro*A) ) ; 


plot (xl , Xnx ( 5 , 

) , 

' k- ' ) 

fn (k) =wn (k) / (2*pi) ; 



plot (xl , Xnx ( 6 , 

) , 

'c-' ) 

fprintf ( ' Mode shape # %2f corresponds 

to nat . 

grid 



freq (fn) : %3.3f\n', k, fn(k) ) 

• 

r 


title ('Mode shapes 

of the Cantilever beam') 

k=k+l ; 



legend ('Mode #1 , 

'Mode #2', 'Mode #3', 'Mode 

end 



#4', 'Mode #5', 

'Mode #6', 0) 


x=linspace ( 0 , L, 180); 
xl=x . /L; 

sigmaN=zeros ( 1 , HMMS) ; 
for ii=l : HMMS ; 


xlabel (' x/L ') ; ylabel('Mode shape X_n(x)') 
holdof f 


end 

end 


sigmaN (ii) = (sinh (betaN (ii) *L) -sin (betaN (ii) *L 
) ) / (cosh (betaN (ii) *L) +cos (betaN (ii) *L) ) ; 

end 

Tc= ' (cosh (betaN (ii) . *x ( j j ) ) -cos (betaN (ii) . *x ( 
j j ) ) ) -sigmaN ( ii ) . * (sinh (betaN (ii) . *x ( j j ) ) -sin 
(betaN (ii) *x ( j j ) ) ) ' ; 

Xnx=zeros (length (betaN) , length (x) ) ; 

for ii=l : length (betaN) 
f or j j =1 : length (x) 

Xnx (ii , j j ) =eval (To) ; 

end 

end 

o 

o 

% Plot mode shapes; 

o 

o 

MMS=HMMS ; 
if MMS==1 

plot (xl , Xnx ( 1 , :) , 'b-') 

title (’Mode shape of the Cantilever beam’) 
legend ('Mode #1', 0); xlabel (' x/L ') ; 
ylabel('Mode shape X_n(x) '); grid 
holdof f 
elseif MMS==2 

plot (xl , Xnx ( 1 ,:) , 'b-'); hold on 

plot (xl , Xnx ( 2 , : ) , 'r-');grid 

title ('Mode shapes of the Cantilever beam') 

legend ('Mode #1', 'Mode #2', 0) 

xlabel (' x/L ') ; ylabel('Mode shape X n(x)') 


II. Two Layer Cantilever Beam 

f unctionCbeam (~) 

% Cbeam.m Cantilever Beam calculations 
% HELP: This script computes mode shapes and 

corresponding natural for 
% beam with two layers (coated beam) 

o 

o 

clearall ; 
clc; 

closeall ; 
display ( ' 

-k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k i 'j 

display ( ' 

~k ~k 'k-k'k'k'k'k'k-k'k-k'k'k'k-k'k'k'k-k'k'k'k'k'k'k > 'j 

display (' Mohammad Abu-Shams ' ) 

display ( ' School fo Engineering and Technology' ) 

display ( ' 

-k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k i ^ 

display ( ' 

-k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k i ^ 

display (' This program computes mode shapes and 
crossponding natural') 

display (' frequencies up to 6 modes for a two 
layers cantilever beam, ' ) 
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display ( ' (the main beam (s) is coated with a thin 
film layer ( f ) ) . ' ) 

display (' Mechanical properties and geometry 
size of the coated cantilever beam') 
display (' should be specified by the user’) 
display (' Prepare to input the following 
proparities : Length, width, thickness (Thf ) , ' ) 
display (' thickness (Ths ), Young Modulus (Es and 
Ef), density (Ros and Rof)') 

display (’The cross section should be 
rectangular, enter 1 if this is correct, other 
wise the program will be closed') 

CS=input (' Enter your choice: '); 

ooooooooooooooooooooooooooooooooooooooooooooo 

oooooooooooooooooooooooooooooooo 

9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 

ooooooooooooooooooooooooooooooooooooooooooooo 

oooooooooooooooooooooooooooooooo 


if CS==1 

W=input ( 'Enter Width of the beam in [m] : ' ) ; 
Ths=input (' Enter Thickness of the main beam (s) 
in [m] : ' ) ; 

Thf =input (' Enter Thickness of the thin film 
layer ( f ) in [m] : ' ) ; 

L=input (' Enter Length of the whole beam in 
[m] : ' ) ; 

As=W*Ths ; 

Af=W*Thf ; 

disp (' Material properties of the beam') 
disp('Do you know your beam"s material 
properties, Young"s modulus and density ?') 

YA=input (' Enter 1 , if you do; enter 0, if you 
dont : ' ) ; 

if YA==1 

Es=input (' Enter Young"s modulus in [Pa] for (s) 
beam : ' ) ; 

Ros=input (' Enter material density in [kg/m A 3] 
for ( s ) beam : '); 

Ef=input (' Enter Young"s modulus in [Pa] for (f) 
beam : '); 

Rof =input (' Enter material density in [kg/m A 3] 

for ( f ) beam : ' ) ; 

else 

display (' Steel : E=2.1e+ll [Pa]; Ro=7850 
[Kg/m A 3 ] ' ) ; 

display (' Copper : E=1.2e+ll [Pa]; Ro=8933 

[Kg/m A 3 ] ' ) ; 

display (' Aluminum : E=0.69e+ll [Pa]; Ro=2700 
[Kg/m A 3 ] ' ) ; 

display (' Zinc-oxide : E=l.lle+ll [Pa]; Ro=5680 
[Kg/m A 3 ] ' ) ; 

display (' PZT : E=0 . 63-1 . le+11 [Pa]; Ro=7750 
[Kg/m A 3 ] ' ) ; 

Es=input (' Enter Young" s modulus in [Pa] for (s) 
beam : ' ) ; 

Ros=input (' Enter material density in [kg/m A 3] 
for ( s ) beam : ' ) ; 

Ef=input (' Enter Young" s modulus in [Pa] for (f) 
beam : '); 

Rof =input (' Enter material density in [kg/m A 3] 

for ( f ) beam : ' ) ; 

end 

else 


F=warndlg ( 'WARNING: INVALID ENTRY !!!','!! 
Warning ! ! ' ) ; 
waitf or ( F) 

display (' Type in:>>Cbeam') 

pause ( 3 ) 

return 

end 


y 

= (Es*Ths A 2+Ef* (2*Ths*Thf+Thf A 2) ) / ( (2*Es*Ths) + 
( 2 *Ef *Thf ) ) ; 

Is = ( (1/12) *W* (Ths A 3) ) + 

W*Ths* (y- (Ths/2) ) A 2; 


If= ( (1/12) *W* (Thf A 3) ) + 
W*Thf* ( (Thf/2) +Ths-y) A 2; 


9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 9 - 

ooooooooooooooooooooooooooooooooooooooooooooo 

OO on ^ / i\/r _m__ 1 'v_'v.0000000000000000000 

o o o o\\l v l(J(ico// ooooooooooooooooooo 

display ( 'How many modes and mode shapes would 
you like to evaluate ? ' ) 

HMMS=input (' Enter the number of modes and mode 
shapes to computed: '); 

if HMMS>=7 
disp ( ' ' ) 

warning (' NOTE : ONLY 6 mode shapes (plots) are 
displayed via the script.'); 
disp ( ' ' ) 

end 

Nm=3*HMMS ; 

j j=i; 

while j j <=Nm; 

betaNL (jj)=fzero(@ (betaNL) cosh (betaNL) *cos (be 
taNL) +1, [ j j j j+3] ) ; 
j j=j j+3; 
end 

index= (betaNL~=0 ) ; 
betaNLall= (betaNL (index) ) ' ; 

%fprintf ( ' betaNL value is %2 . 3f \n ' , betaNLall ) ; 
betaN= (betaNLall/L) '; 
k=l ; 

wn=ones (1, length (betaN) ) ; 
fn=ones (1, length (wn) ) ; 
while k<=length (betaN) ; 

wn (k) =betaN (k) A 2* (sqrt ( (Es*Is+Ef*If) / (Ros*As+ 

Rof*Af ) ) ) ; 

fn (k) =wn (k) / (2*pi) ; 

fprintf ( ' Mode shape # %2f corresponds to nat . 
freq (fn): %3.3f\n', k, fn(k) ); 
k=k+l ; 

end 

x=linspace ( 0 , L, 180); 
xl=x . /L; 

sigmaN=zeros ( 1 , HMMS) ; 
for ii=l : HMMS ; 

sigmaN (ii) = (sinh ( (betaN (ii) *L) ) -sin (betaN (ii) 
*L) ) / (cosh (betaN (ii) *L) +cos (betaN (ii) *L) ) ; 

end 

Tc=' (cosh(betaN(ii) ,*x(jj))-cos(betaN(ii) . *x ( 
j j ) ) ) -sigmaN (ii) . * (sinh (betaN ( ii ) . *x ( j j ) ) -sin 
(betaN (ii) *x ( j j ) ) ) ' ; 

Xnx=zeros (length (betaN) , length (x) ) ; 

for ii=l : length (betaN) 
f or j j =1 : length (x) 

Xnx ( ii , j j ) =eval (Tc) ; 

end 

end 

o 

o 

%Plot mode shapes; 

o 

o 

MMS=HMMS ; 
if MMS==1 

plot (xl , Xnx ( 1 ,:) , 'b-') 

title ('Mode shape of the coated Cantilever 
beam ' ) 

legend ('Mode #1', 0); xlabel ( ' x/L ' ) ; 
ylabel('Mode shape X_n(x) '); grid 
holdof f 
elseif MMS==2 

plot (xl , Xnx ( 1 ,:) , 'b-'); hold on 

plot (xl , Xnx ( 2 , : ) , 'r-');grid 

title ('Mode shapes of the coated Cantilever 

beam ' ) 

legend ('Mode #1', 'Mode #2', 0) 

xlabel (' x/L ') ; ylabel('Mode shape X_n(x)') 

holdof f; 

elseif MMS==3 
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plot (xl , Xnx ( 1 ,:) , 'b-'); hold on 

plot (xl, Xnx (2 , : ) , ' r- ' ) 

plot (xl , Xnx ( 3 , : ) , 'm-');grid 

title (’Mode shapes of the coated Cantilever 

beam ' ) 

legend ('Mode #1', 'Mode #2', 'Mode #3', 0) 
xlabel ( ' x/L ' ) ; ylabel('Mode shape X_n(x)') 
hoi doff; 
elseif MMS==4 

plot (xl , Xnx ( 1 , :) , 'b-'); hold on 

plot (xl, Xnx (2 , : ) , ' r- ' ) 

plot (xl , Xnx ( 3 , : ) , 'm-') 

plot (xl , Xnx ( 4 , : ) , ' c- ' ) ; grid 

title ('Mode shapes of the coated Cantilever 

beam ' ) 

legend ('Mode #1', 'Mode #2', 'Mode #3', 'Mode 
#4', 0) 

xlabel (' x/L ') ; ylabel('Mode shape X_n(x)') 
hoi doff; 
elseif MMS==5 

plot (xl , Xnx ( 1 , 'b-'); hold on 

plot (xl, Xnx (2 , : ) , ' r- ' ) 

plot (xl , Xnx ( 3 , : ) , 'm-') 

plot (xl , Xnx ( 4 , : ) , ' g- ' ) 

plot (xl , Xnx ( 5 , : ) , ' k- ' ) 

grid 

title ('Mode shapes of the coated Cantilever 
beam ' ) 

legend ('Mode #1', 'Mode #2', 'Mode #3', 'Mode 
#4 ' , 'Mode #5 ' , 0) 

xlabel (' x/L ') ; ylabel('Mode shape X_n(x)') 
holdof f 
elseif MMS>=6 

plot (xl , Xnx ( 1 ,:) , 'b-'); hold on 

plot (xl, Xnx (2 , : ) , ' r- ' ) 

plot (xl , Xnx ( 3 , : ) , 'm-') 

plot (xl , Xnx ( 4 , : ) , ' g- ' ) 

plot (xl , Xnx ( 5 , : ) , ' k- ' ) 

plot (xl , Xnx ( 6 , : ) , ' c- ' ) 

grid 

title ('Mode shapes of the coated Cantilever 
beam ' ) 

legend ('Mode #1', 'Mode #2', 'Mode #3', 'Mode 
#4', 'Mode #5', 'Mode #6', 0) 
xlabel (' x/L ') ; ylabel('Mode shape X_n(x)') 
holdof f 


end 


end 
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